Nonlinear waves of nuclear density 



V.G. Kartavenko a ' c ' 1 , A. Sandulescu b and W. Greiner c 

a Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear 
Research, Dubna, Moscow District, 141980, Russia 

^Romanian Academy, Calea Victoriei 125, Bucharest, 71102, Romania 

c Institut fur Theoretische Physik der J. W. Goethe Universitat 
D-60054 Frankfurt am Main, Germany 



Abstract 

Nonlinear excitations of nuclear density are considered in the framework of semiclas- 
sical nonlinear nuclear hydrodynamics. Possible types of stationary nonlinear waves 
in nuclear media are analysed using Nonlinear Schrodinger equation of fifth order 
and classified using a simple mechanical picture. It is shown that a rich spectrum 
of nonlinear oscillations in one-dimensional nuclear medium exist. 
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1 Introduction 



Nonlinear aspects of clusterization, one of the most mysterius and important 
problem of modern nuclear physics, are of interest due to the following reasons. 

First, practically all nuclear processes which are related to clustering lead 
to large reconstruction of an initial nuclear system and this requires a de- 
velopment of new theoretical methods to describe large amplitude collective 
motions. 

Second, the existance of clusters is a very general phenomenon. There are 
cluster objects in subnuclear and macro physics. Very different theoretical 
methods were developed in these fields. However, there are only few basic 
physical ideas, and most of methods deal with nonlinear partial differential 
equations. 

Third, methods of nonlinear dynamics give us possibility to derive for nuclear 
physics unexpected collective modes, which can not be obtained by traditional 
methods of perturbation theory near some equilibrium state. 

The alpha decay is the most studied process among the other possible channels 
of fragmentations at low energies, such as cluster radioactivity, cold fission and 
fusion[l]. As a rule these processes lead to a large mass transfer up to the final 
channel. The existence of a central region of a constant density and a well 
shaped surface region makes it possible to describe a clusterization in the 
region of nuclear surface including the neck region [2]. 

The alpha decay and cluster radioactivity indicate well the formation of a clus- 
ter on the surface of a large nucleus [1]. Our nonlinear approach was formulated 
in such a way to describe the preformation of such cluster [3]. Our model is 
an extension of the Bohr-Mottelson collective model[4], which allows not only 
the nuclear deformations which lead to collective vibrational and rotational 
motions, but also the creation of bumps on nuclear surface. 

Nuclear density falls considerably down in a region of nuclear surface. This fact 
is very important due to possibility of a clusterization and other fluctuations 
of a nuclear density leading to instability in a surface or a neck regions [3], [5]. 
Recently[6] we suggested an alternative nonlinear approach, based on the cur- 
rents and density algebra [7] to describe a surface interaction of a nuclei with a 
cluster. Starting from the hydrodinamic set of equations and using an effective 
Skyrme interaction, a nonlinear Schrodinger equation with fifth order term in 
the nonlinearity (NoSE3 — 5) has been deduced to describe in a semiclassical 
limit an irrotational flux of a nuclear density. In the framework of this ap- 
proach an approximate way to describe surface axial symmetric interactions 
of a large target with a small cluster has been derived. 
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Fig. 1. Surface interaction of 08 Pb and Mg, as a superposition of nonlinear waves 
of NoSE3-5 

Figure 1 gives the density distributions for a surface axial symmetric (along 
the z-ax.es) interaction of a large target with a small cluster, as an interference 
of two nonlinear solutions of NoSE3-5. One can see the transition from well 
localized and separated initial distribution to a complicated picture, when the 
nucleus 28 Mg begins to be absorbed in the surface region. 

Methods of nonlinear dynamics are rather unusual for traditional nuclear the- 
ory. Therefore it would be useful to understand: Why there exist nonlinear 
effects within linear Schrodinger picture with a many-body nuclear hamilto- 
nian? Does there exist a way to obtain nonlinear evolution equations, starting 
with traditional Hamilton or Lagrange picture? How to classify possible types 
of solutions and to understand their physical meaning? 

Within this paper we shall try to clarify these points aimed to consider simplest 
large amplitude oscillations of nuclear density. 

The paper first presents in Sec. 2 the basic scheme of our nonlinear approach 
to nuclear hydrodynamics. In Sec. 3 it is shown how to reduce a set of hy- 
drodynamical equations to a dimensionless nonlinear Schrodinger equation. 
Properties of nonlinear stationary waves are investigated in Sec. 4. We give the 
answers to these three above questions and present a short summary in the 
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last section 



2 Nonlinear approach to nuclear hydrodynamics 

Fluid dynamical concept [8] of the nuclear medium is a traditional way to clas- 
sify collective modes. Hydrodymamics gives up natural set of variables(collective 
currents and density) to describe any collective excitations of nuclear density. 
At present there are a lot of formulations of the nuclear hydrodynamical ap- 
proach. We will follow the way, which is based on the current and density 
operators algebra [2], [7] 

2. 1 General scheme 

For nuclear systems a usual second-quantized form of a nonrelativistic theory is 
defined by introducing canonically conjugated nucleon fields ^(x, a, q), where 
(q) denotes (n,p) for neutrons and protons and a — ± the spin index, which 
satisfy the equal-time canonical anticommutation relations 

{* + (l), *(2)} + = 5(1 - 2), {*+(!), * + (2)} + = *(2)} + = 0.(1) 

In terms of these operators a nonrelativistic Hamiltonian could be defined 

H nucl [^ + ^- U nud {x)\ (2) 

with a general two-body nuclear interaction U nuc i(x) 

The term "hydrodynamics" means that we will describe the dynamical behav- 
ior of the nuclear system in a restricted space of collective variables represent- 
ing the density and the nucleon current of the system: 

(*+(£), =► },(£)) ( 3 ) 

where the index k denotes the cartesian vector components. 

The description of the evolution in terms of collective currents (3) requires 
the of use commutation relations between the collective current and density 
operators 

[p(x),jk(y)}, [p{x),p{y)}, \jk(x),]i(y)}, (4) 
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and their commutation relations with the hamiltonian (2) 

Op(xj 1 /v djfc(x) 1 

= \p(x) , H nud ] , — = — \j k (x) , H nud ] . (5) 



Commutation relations (4,5) can be obtained by using the initial commutation 
relations (1) for nucleon quantum fields and the existing formal representation 
for the kinetic energy density tensor operator 

f nk (x) = V n V + (x) ■ V k V(x) + V k V + (x) • V n *(£) (6) 

in terms of collective currents 

TT = (hp-'ji + JiP^Jk) + \p~ l V fe p • V lP + const • 5 kl (7) 

which satisfies Eqs. (5) as well. 

As a result one gets a collective hydrodynamical hamiltonian 

H nud [V + ; *; U nucl (x)] J; T^, <S[p, /]] (8) 

which is equivalent to the initial nuclear hamiltonian (2) relative to the equa- 
tions of motion (5). Here <B[p, j] denotes effective potential term. We present 
the hamiltonian only schematically aimed to show the general scheme. 

It is convenient to use Skyrme interaction, which gives an ansatz for an effec- 
tive hamiltonian Hhyd / d 3 xH(x) of a local operator in terms of densities 
and currents [1 



H(x) = (h 2 /2m)r + B lP 2 + B 2 {pl + # + B 3 (pr - f ) 

+ B 4 (p n r n - 2 + ppT p -T P ) + B 5 pAp + B e (p n Ap n + ppAp p ) 

+ B 7 p 2+a + B 8 p a (p 2 n + p 2 p ) + S 9 (pV • S + J • V x E 

+ p n V • S„ + j n • V x S n + p P V • S p + j p • V x S p ) 

+ B w f? + B n (S + Ej) + 5 12 p a E 2 + 5 13 p a (^ + gj), (9) 

where the numerical coefficients B t and a are parameters of a Skyrme inter- 
action [9]. 

To describe a system of nonrelativistic spin-1/2 nucleons in terms of currents 
requires to decompose spin-isospin density tensors to spin-scalar 

p,(f)EE^H/ + (f,a,g)vl/(f,a,g) (10) 
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and spin-vector densities 

= J2 a > s) < °\°W > > ^ ?)■ ( n ) 

<J,a' 

The related currents are the spin-scalar 

Jk(x) q = ^ ^ ?)V fc ^(x, a, g) - V fc M/ + (f , a, q) • a, g)) , (12) 

and spin-orbital ones 



Sfc(x) g = S ^kmn(^ + (x,(r,q)V m ^(x,a',q) (13) 

ct,ct' m,n=l 

- V m y + (x,a,q)V(x,a',q)) < a\a n \a' > . 

Among all components of the kinetic-energy-spin-tensor the most important 
ones are the scalars in spin and in the ordinary coordinate space 

r q , j 2 q , V-S„ (14) 

1 ~ _ ~ 1 

T = Yl T kk = S ifcP" 1 ^ + TP" 1 VfepV fe p. (15) 
Z fc fc 4 

It should be noted that the operators (14-15) is not the complete set of spin- 
isospin tensor components of the current-spin-tensor and kinetic-energy-spin- 
tensor. The detail analysis of the equation of motion (5) for all defined densities 
and currents could help to clarify the situation. First of all it is necessary to 
define the effective hamiltonian. 

In order to analyze the evolution equations (5) with the effective hydrody- 
namical hamiltonian (9) it is necessary: first to derive commutation relations 
between all components of spin-isospin current and density tensors, and their 
commutation relations of with all terms of hamiltonian (9) and second to pro- 
vide tensor decomposition of all currents to irreducible cartesian components 
aimed to separate irrotational and rotational terms. 

This task is rather difficult. However, in the framework of presented above 
general quantum hydrodynamics, one could build, in principle, a quantum 
approach to investigate spin-scalar, spin-vector, spin-orbital vibrational and 
rotational modes in nuclei. Such investigations are in progress. 
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3 Nonlinear Schrodinger equation 



In order to check the possibilities of our approach we simplify the hamiltonian 
(9) aimed to describe different oscillations of isoscalar nuclear density. We 
assume a system of spinless and isospinless nucleons with no Coulomb, spin 
and spin-orbital effects (W = 0) and effective Skyrme forces (9) (a = 1) 



H W = ^ r + ^op 2 + ^3P 3 +) + ^(9tx - 5t 2 )|Vpf 



+ -(M 1 + 5t 2 )({p,r} + -2j 



7*2 



(16) 



Schematic " nonlinear density and current dependence" of kinetic energy is 

TT _ ft 2 3 .1. frVlVpl 2 3 2 1 , 

2m* j—t p 8m* p 8 16 



where the following renormalization is used 



m ■■ 
8m 



m = 
H 2 r] 2 
8m 



{m- 1 + (3ti + 5t 2 )p 7V /8^ i ) 



2\-l 



8m 64 v 



5t 2 



-f^mod form Eq. (17) is the simplest hamiltonian which includes the "collec- 
tive" kinetic energy, "Weizsacker" kinetic part which includes the quantum 
and other gradient terms and the "potential" energy, containing effective two- 
body attractive term and three-body imitation of Pauli principle. 

In classical hydrodynamics [10], the fundamental variables are the local density 
p(x) and the fluid velocity v(x) from which the quantum current density can 
be defined by the anticommutator[ll][^| 

in 

= ^{p(x),Vi(x)} + . (19) 
Using the definition (19) and the commutative relations between density and 



2 See discussion on a velocity concept in quantum hydrodynamics in a recent pa- 
per [14] 
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current operators (4) one derives velocity-density and velocity-velocity com- 
mutative relations 

[vk(x),p(y]}, [p(x),p(y]}, (20) 

To classify the type of motions in the cartesian coordinate space it would be 
convenient to provide tensor decomposition of the operator vector field 

v k (x) = V k % ot (x) + curl k (A), div(A) = 0, (21) 

into a sum of potential (& pot (x)) and solenoidal (A) fields by using the Helmholtz 
theorem (see, e.g., ref. [12] vol. 1). 

In the following we shall consider the equations of motion (5) with the hamil- 
tonian H mod (17). First it is convenient to use the scale transformations [7] 

x = (hr]/m*c s )x\ t= (m*c 2 Jh)T, (22) 

v k (x) = c s U k {x') , p(x , t) = p N n 2 (x 1 , t) . (23) 

where c s is the sound velocity in nuclear matter 

= 1 2 5\<E/p) 1/2 
s - [ m* P 5pi 1p=Pn) 

The scale transformation (22,23) makes it possible to separate the general 
effects related to nonlinearity in the hamiltonian (17) (a polinomial in density) 
from effects related to a choice of parameters of interaction, which define only 
scales factors: effective length wave (Hr)/m*c s ), the density of nuclear matter 
Pat, the sound velocity in nuclear matter c s and the time factor {m*c 2 s /h). 

Hamiltonian (17) can be cast in the following form 

Now the continuity and Euler operator equations (5) can be reduced to the 
following set of the dimensionless equations of motion 
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(24) 



hj-± = - (([/ x cw/f7) fe - (curlU x [7) J 

7 -si^^'-^-^+H- (27) 



Taking into account the velocity-velocity and velocity-density commutative 
relations (19-20) and the operator Helmholtz theorem (21) we obtain the fol- 
lowing decomposition of the dimensionless velocity operator: 

d®(x) - ~ ~ 

U k (x) = — ^ + R k (x), divR = 0, C = curlR, (28) 

OXh 



into potential and pure rotational currents. 



Considering only the vibrational type of motion and neglecting their possible 

connection with rotational modes (( = 0, Rk — > 0) the commutation relations 
between density and potential velocity are reduced to the canonical boson 
form 

[n(x)My)]=M@-y), = 0, [h(x),h(y)} = 0, (29) 



which give it possible to treat (if it is no rotation (R k = 0)) the evolution via 
Hmod as nonlinear operator anharmonical picture. 

In the semiclassical limit the two hydrodynamical equations of motion (26,27) 
in the case of irrotational flow (28,29) can be cast in the form of one nonlinear 
Schrodinger equation 

iip T = -Atp - 4|^| 2 ^ + 3|^| 4 ^ (30) 



for a complex function ip 

ti(x,t) =\ip(x,r) | 2 , $(x, r) — argip(x, r). (31) 



Eqs. (30,31) become transparent if we recall the well-known formal analogy of 
quantum mechanics to fluid mechanics[8]. However one should keep in mind 
that in formulaes (30,31) all variables are collective ones and the effective 
potential is defined by a first functional derivative of effective interaction on 
density. 
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4 The one-dimensional nonlinear waves 



Let us consider the stationary (travelling) solutions of the one-dimension (in 
the cartesian coordinate space) nonlinear Schrodinger equation (30). These 
solutions take the usual form 

ip(r, r) = ip(x) exp(i(A — v 2 /2)r + ixv + cr ) (32) 



where r = x — vr — r . The parameters r , ctq, v are the initial position of the 
center of the wave, the initial phase shift and the initial velocity. 

The function ip(x) is the solution of the nonlinear Schrodinger equation 

i) xx -\i/> + 4^ 3 - 3ij 5 = 0. (33) 




Fig. 2. Dependence of "potential" V(ip) on ip and and A 



To analyze different solutions of Eq. (33) we will exploit its simple mechanical 
analog, in which a macroscopic point "particle" of unit mass is moving under 
a conservative "force" 

F(V>) — Xip- 4^ 3 + 3^ 5 = (34) 
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where the "potential" (see Figure 2) is taken as 



(35) 



with if; as a "particle position" and x its "time". 

Eq. (33) after being multipled by ip x , can be one time integrated. The result 
is 

i> x = ±y/2(E-Vty)) t E = + X -^ x (36) 

where the constant of integration (E) is considered to be the "energy" . Eq. (36) 
is just the energy conservation law in the analog "particle" problem. The 
parameter (A) defines the "potential" profile and the "energy" ( E ) defines 
the "particle" motion. 

Figure 2 gives a gross "potential" picture. We show a detail set of "potentials" 
V(ip) for different values of the parameter (A ) in Figs. 3 . The "potential" 
V(ip) (35) vanishes at the points 

i) = {0, ±Pi = ±y/l- VT^J, ±P 2 = ±y/l + VT^J} (37) 
(see Figure 3 (A = 0.5)). 

The equilibrium points of motion F{iji) = —dV{jp)/dip = are given by 

, ={0 , ±P3 = ±f-^qEE, ±Pi=± ^±^ } (38) 

(see Figs.3 (A = 0.9)). For an nonpositive A < a finite amplitude motion is 
possible only for < E < V(P 4 ) and |^(0)| < P 4 (see Figs.3 (A = 0.0)). For 
< A < 1 there exist a finite amplitude motion eith different types for the 
three regions of the "energy" {V(P 3 ) < E < 0, E = 0, < E < V(P 4 )} (see 
Figs.3 (A = 0.5, 0.9)). This is the most interesting case and in the following we 
will consider it in detail. For A = 1 the "barriers" V(-P A ) = V(0) = V(P A ) = 
(Figs.3 (A = 1.0)). A finite amplitude motion for 1 < A < 4/3 is possible only 
for the negative "energy" V(P 3 ) < E < (Figs.3 (A = 1.1)). For A > 4/3 (see 
Figs. 3 (A = 1.5)) any motion of a finite amplitude is impossible. A "particle" 
comes away from the origin for any energy. 

Eq. (33) is an ordinary differential equation of a second order. For any set of 
initial values ^(^o); i>x( x o) & t initial "time" x one can obtain, in principle, the 
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solution of Eq. (33) numerically. We will restrict ourselves only to stationary 
waves of a finite amplitude. This means that for a given A the "energy" is less 
than the maximum "barrier" value ( E < V(±P^)). Therefore it is convenient 
to select an "initial point" x , as the point at "potential" profile in which 

E = V{ip ), ^(x ) = 0, ij = ij(x ). (39) 



Because of the "energy conservation" the E parameter is related to the wave 
amplitude ip(0). A localized soliton like solution is the nodeless solution of 
Eq. (33) under the following boundary condition 

V>(oo) =^(oo) =^(0) = 0. (40) 



This corresponds to a special case with E = 0. Soliton like solution exists if 
< A < 1. A "particle" starts motion in the point P\ or —Pi (see Eq. (37)) 
and comes to origin. 

In Figures 4 we present the soliton like solutions for a few values of the pa- 
rameter (A). There exist an analytical form of this soliton like solution[7] 

rl){x) = (A/(l + (1 - A) 1/2 cosh(2A 1/2 r)) 1/2 . (41) 



which looks like a symmetrized Fermi-distribution 

ip(x) 



sinh(^/D) 



where the "radius" 



\ cosh(9fy2)) + cosh(x/D) 



V\= ^ 7 =arctanh(v / A). (43) 



and the " diffuseness" 



s = irx (44) 



are defined by the parameter (A). 

The parameter A is the chemical potential defined by the "particle number" 

oo 

N= J dxip(x) 2 , N = arctanh(VT) (45) 

— oo 
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Fig. 3. Dependence of a "potential" V(tp) on ip for different parameters A 
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Fig. 4. Soliton -like solutions of NoSE3-5 for different values of the parameters A 

It is useful to note also the following features of this soliton like solution. 
First for a small "particle number" the solution looks like a solution of a 
cubic nonlinear Schrodinger equation. Second, for an intermediate and large 
"particle numbers" the "radius" (43) is increased (see Figs.4). Third, there 
exist an internal region of approximate constant "density" and a "surface 
region" of a "constant" diffuseness (see Eqs. (43,44) and Figs.4(c,d)). 

In Figures 5 we present possible periodical and aperiodical solutions of Eq. (33) 
for A = 0.5. The corresponding "potential" is represented in Figs. 3 (A = 
0.5). One can see small oscillations with a positive amplitude in the "well" 
for (V(P 3 ) < E < 0) (Figs. 5(a)) These oscillations became distorted more 
and more, as E — > — e, where e is an infinite small positive value (see 
Figs.5(b,c)). For E = we have the exact soliton like solution (Figs. 5(d)). 
There are nonlinear waves of positive and negative apmlitude, as E — > V(P 3 ) + 
e, (see Figs.5(e,f)). These waves look like quasi-trigonometrical functions (see 
Figs.5(g)), as A changes from small positive values to the "barrier" V{P±). 

Periodical waves became distorted more and more, as E — > V(P±) — e. In 
Figs.5(h,i) one can see a transition from a well-defined periodical structure to 
an aperiodical one. 
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Fig. 5. Set of stationary nonlinear waves of NoSE3-5 
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There is a region of semi "shock" wave with the well defined sharp "front" (see 
Figs.5(i)). "Front" moves to the infinity and, in this case, we have the constant 
solution in coordinate space when an "energy" is equal to the "barrier" (E = 
V(P 4 )) (see Fig.5(j)), describes a one-dimensional nuclear matter. 

The above solutions were obtained numerically. However it is possible to get 
also analytical solutions for few special cases. First it shoud be noted that it 
is possible to write a formal solution of Eq. (36) 

ip dip , 
± xdx 



J xdx (46) 



for any types of nonlinearity, i.e. for any "potential" V{ip). 

In our case V(ip) is given by a symmetrical polinomial of six order in ip (35). 
The substitution Z = ip 2 reduces equation (33) to the following form 



Z x = ±2^Jz(2E + XZ - 2Z 2 + Z 3 ) 



= ±2^Z(Z - Z l )(Z - Z 2 )(Z - Z 3 ) (47) 

The expression under the square root of the Eq. (47) is a polinomial of fourth 
degree. One of the roots is Z = the other three ones being defined by E, A 
parameters. 

The following solution of Eq. (47) exists also in the case of E = 0, A = 

JZ{Z - 2) 2 

= ±2x Z(x) = (48) 

2Z y ! l-Ax 2 K J 



where the constants of integration are omitted. This solution has a pole at 
x = ±2, and it does not describe a finite solution. 

For very few special cases Eq. (47) can be integrated analytically in Jaco- 
bian elliptic functions. These functions are quite similar to the trigonometric 
functions. One can find in the book[13] the complete set of twelve different Ja- 
cobian ellipic functions and the coresponding form of the differential equation. 
If it is possible to reduce Eq. (47) by a linear transformation (Z, x) — * (z, u) 
to the one of these forms, then the solution of (47) can be cast as the related 
Jacobi elliptic function. The explicit form of the solution is determined by the 
constants A, E and boundary conditions. 
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5 Summary 



We give the following answers to the three questions formulated in Introduc- 
tion. 

Nonlinear problems exist within a linear Schrodinger picture with a many- 
body nuclear hamiltonian as the result of the reduced description in a re- 
stricted space of collective degrees of freedom. In Sections 2 and 3 the way 
to derive the Nonlinear Schrodinger equation for a collective motion is pre- 
sented. The logic of our approach corresponds to the logic of any microscopic 
or semimicroscopic theory of collective nuclear motion, namely, one first selects 
a space of certain collective variables (3) and their commutation relations (4). 
Then seeks an expresssion for the collective Hamiltonian (8) in terms of the 
chosen collective variables to reproduce the commutation relations (5) of the 
original many-particle Hamiltonian and the collective operators. The hydro- 
dynamical collective hamiltonian (8) is equivalent to the initial Hamiltonian 
(2) only with respect to the equations of motion (5) for the collective operators 
(2). The present approach is the realization for nuclear physics of the general 
problem how to formulate nonrelativistic quantum mechanics in terms of col- 
lective current operators. The currents (3) can be considered as fundamental 
dynamical variables rather than the underlying canonical fields (1). The main 
emphasis of this type of theory is the algebraic structure of the equal-time 
current-current commutation relations (4). The motivation for initiating such 
a program was due to the fact that most of the interesting physical quantities 
can be expressed in terms of currents and densities. It is important for nuclear 
physics because the density and velocity distributions are usual variables to 
analyse the dynamics of heavy ion collisions (see the recent papers [14], [15] and 
the references therein). 

In the present paper, possible types of stationary nonlinear waves in nuclear 
media are analyzed in the framework of Nonlinear Schrodinger equation of 
fifth order and classified using a simple mechanical picture. It is shown that 
a rich spectrum of nonlinear oscillations in one-dimensional isoscalar nuclear 
matter exist. 
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